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Abstract: If R-parity conserving super symmetry is realised with masses below the TeV 
scale, sparticles will be produced and decay in cascades at the LHC. In the case of a neutral 
LSP, which will not be detected, decay chains cannot be fully reconstructed, complicating 
the mass determination of the new particles. In this paper we extend the method of 
obtaining masses from kinematical endpoints to include a gluino at the head of a five- 
sparticle decay chain. This represents a non-trivial extension of the corresponding method 
for the squark decay chain. We calculate the endpoints of the new distributions and assess 
their applicability by examining the theoretical distributions for a variety of mass scenarios. 
The precision with which the gluino mass can be determined by this method is investigated 
for the mSUGRA point SPS la. Finally we estimate the improvement obtained from adding 
a Linear Collider measurement of the LSP mass. 
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1. Introduction 



Supersymmetry p], |2l ||, || provides one of the more popular extensions to the Standard 
Model at higher energies. It has many attractive features, one of which is a possible solution 
to the hierarchy problem ||. For this to be the case supersymmetric particles must exist 
near the TeV-scale, and may therefore be accessed by the LHC. 

Conservation of R-parity means that any interaction vertex must involve an even num- 
ber of supersymmetric fields. For collider experiments this has two important consequences: 
sparticles will be produced in pairs, and the decay chain of a sparticle will always end with 
the lightest supersymmetric particle (LSP). A SUSY event at the LHC will then result in 
two decay chains, each giving an undetected LSP together with a number of detectable 
Standard Model particles. The escaping LSPs make it difficult to fully reconstruct events. 
This means that the masses of the sparticles involved in a decay chain cannot be readily 
reconstructed from the end products. However, even though they do not represent any 
particular sparticle, the various mass distributions that can be constructed from the de- 
tectable particles of the decay do have a sensitivity to the sparticle masses. In particular, 
the maximum value of each of these distributions can in principle be calculated for any 
given decay chain, and this maximum value will be determined by the masses involved 
in the decay. If these kinematical endpoints are measured, it is then possible to obtain 
the sparticle masses in a numerical fit. This 'endpoint method' has been widely used to 
determine masses of supersymmetric particles || [7|, ||, |l0|, 11 ] . The decay chain 



is particularly amenable to this method 1 . In [12| we investigated this decay chain in 



considerable detail for two points on the SPS la line |13|]: (a), which is the standard point 
and often simply referred to as 'the SPS la point', and which lies at somewhat higher 
masses. In this paper we restrict ourselves to the standard SPS la point and investigate 
the situation where a gluino is at the head of the decay chain, 

9 -> <7L<?n X2<Mn -> Wn<?f<?n ~> XlVn<Mn (1-2) 

Seven more distributions become available with the inclusion of q n , and their endpoints, if 
measurable, will enable us to find the gluino mass. For ease of reference the new distribu- 
tions/endpoints will be called 'gluino distributions/endpoints' as opposed to the 'squark 
distributions/endpoints' which do not involve the gluino mass. 

Although very much in line with the now standard method for the masses involved 
in (p. . 1|) , this way of obtaining the gluino mass is new. At present two other methods are 



proposed. One is an approximate technique, used for example in [14]. If, in the three-body 



decay ~^ Xv-h we have mu at its maximum value, ran = "i™ ax , then X\ 1S a t rest 



1 Maintaining the standard convention, the subscripts 'n' for near and 'f for far on the leptons and 
quarks are used to distinguish the order of particle emission in the decay chain. These subscripts will be 
suppressed if there is no ambiguity. In |l2| , which will be frequently referred to throughout the text, no 
subscript is used for the quark as there is only one. In the current notation this quark corresponds to qt . 
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in the rest frame of x\- Assuming m^o known, the four- vector of X2 can then be fully 
reconstructed. Next, adding the two quark four- momenta provides us with the complete 
four-momentum of the gluino (and squark). This allows the mass of the gluino (and the 
squark) to be found in the more traditional way of plotting the mass peak. Apart from the 
normal experimental uncertainties this method is exact. However, for the decay chain (|1.2|), 
it is generally not the case that Xi is at rest in the rest frame of X2 when mu = m™ ax . Only 
if m~ = m^om^o is this true, and the more this relation is violated, the less reliable the 
results will become. Furthermore, to get a sizable sample, also events with mu somewhat 
below m™ ax must be used, which increases the uncertainty of the resulting masses. So, to 
apply this method to ( |1.2|) the following three effects must be controlled: i) the knowledge 



of m^o — this can be assumed known from the standard endpoint analysis of the squark 



chain, but the appropriate uncertainties must be included, ii) the systematics from the 
violation of m~ = m^m^o, and iii) the systematics from the inclusion of events with 
mu < m™ ax . A rigorous treatment of all these effects is still wanting. 

The other method to obtain the gluino mass from (1.2) goes under the name of the 
'Mass relation method', and is described in [15|. In a first version the masses of Xi, Ir and 
X% are assumed known, e.g. from the squark-endpoint analysis. For a given gluino-chain 
event we then have six unknowns: mg, ma* and the four- momentum of x±- On-shell mass 
conditions for the five sparticles provide five equations. For one such event we are one 
condition short of determining the system. [If also m^ L is assumed known, the system is 
solvable.] A second event of the same type will again give 5 equations for 6 unknowns. Two 
of the unknowns, m§ and mg L are however the same for the two events. In combination 
one therefore has 10 equations for 10 unknowns. The system can be solved analytically to 
provide (in some cases multiple) masses for the gluino and the squark. The mass relation 
method involves no approximations of the type described for the first method, and is in 
this respect at the same level as the endpoint method we will use here. The mass relation 
method also has the advantage that it avoids the difficulties related to estimating endpoints, 
and it can give a measurement even for quite a small number of events. On the other hand 
the uncertainties on the three lighter masses as well as possible correlations between the 
masses need to be appropriately handled. An extension of the method to take no input 
masses, and/or include in a consistent way measurements from an endpoint analysis should 
be feasible and would meet these requirements. 

For the endpoint method, which in the following is developed to also provide the gluino 
mass, all measurements go into the same fit, so errors and correlations are automatically 
treated correctly. 

The structure of the present paper is as follows. In Sect. [2] we present the endpoints 
of the invariant distributions associated with gluino decays. We also discuss the ambiguity 
which arises from having two quarks in the final state and plot some examples of gluino 
distributions with no experimental effects. Sect. || presents a simulation of the gluino 
endpoints at the ATLAS detector for the benchmark SUSY scenario SPS la, ending with 
an estimate for some of the new endpoints. Energy scale errors are also discussed in some 
detail. In Sect. |] values for masses and mass differences are found from the estimates 
of the endpoint measurements. In Sect. |B| the impact of the statistical error is separated 
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from that of the energy scale error and two somewhat more conservative error scenarios 
are studied. The effect of adding a Linear Collider measurement is investigated in Sect. |6|, 
while conclusions are made in Sect. ^. A new method for calculating rn^f is developed and 
generalised to the gluino case m^f in App. [A|. Finally, the complicated gluino endpoint 
m gg/fiow) 1S derived in some detail in App. [B]. 

2. Gluino endpoints 
2.1 Endpoint calculation 

Using the method detailed in App. [A] we obtain the maximum value of m qq u for the five- 
sparticle gluino chain: 
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Notice the systematic form of the defining inequalities and the resemblence to rn^f^, given 



in Eq. (4.4) of [12]. For the gluino endpoint we have four 'dominance regions', each defined 
in terms of a nearest-neighbour mass ratio dominating the product of the three others. 
As before, simplifications can be made to these inequalities, but with the undesired conse- 
quence of obscuring the neat and systematic structure. 

In principle, a measurement of rn^f would, in combination with the squark endpoints, 
suffice to provide the gluino mass. The long decay chain does however allow for more 
distributions to be constructed, from which endpoints can be extracted and compared with 
analytic expressions. These other available endpoints are useful both as over-constraining 
measurements and as consistency tests. 

The gluino chain has seven primary endpoints which involve nig. The two-particle 
primary endpoints easily calculated from cascade decays 'on a line'. 

There is only one realisation for each. 

For the three-particle endpoint rn™^, the particles involved are neigbours in the decay 
chain, giving a situation similar to w-™^ x . The solution can be found from Eq. (4.4) of |l2| 
via the substitutions 



(m u , m^o , m~ lR , m^o ) -> (mg , m u , m^o , ) 
in both the endpoint expressions and the inequalities. 



(2.2) 
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The two other three-particle endpoints and are more difficult to find be- 

cause the particles involved are not all nearest neighbours in the decay chain. However, 
it is possible to transform the problem into one which involves only nearest neighbours. 
Let us first consider Tn q jf. The initial point of difficulty is the orientation of the decay, 
where the unused qf is emitted. Will the orientation which gives the largest m qn a depend 
on the mass scenario? To answer this question, boost to the rest frame of x% Whatever 
the details of the two steps prior to the creation of jfy, at this stage the only quantity 
of relevance to Tn qn ii is |p g J- The larger the momentum is, the larger m qn i\ can become. 
Thus, independent of masses, for maximum values of m qn ii, the orientation of the q~L decay 
is always the one that maximises [p> | in the rest frame of x\- This means sending x\ 
in the opposite direction of q n , and the configuration which gives rn^f therefore has p qn 
parallel to p,j f . 

With this point settled, the situation can be transformed into a nearest-neighbour 
decay, for which we know the solution. To do this we introduce a pseudo-particle: as seen 
from the rest frame of \\ the momentum of q n corresponds to the decay of a pseudo-particle 
Y — > XzQn with mass my = ni'- — m| £ + m| . In such a picture all particles are nearest 



neighbours, and the solution Eq. (4.4) of p2| applies with the appropriate substitution 



m qL ^ my = ^/m? - m? L + m|o (2-3) 

For the orientation of the X2 decay is fixed to maximise |pj f | in the rest frame 

of x[>! analogous to the argument above. Another pseudo-particle is then defined at the 
end of the decay chain, X2 ~^ wnere m Y = m x\ ra x\l m lR' Again, the problem has 
been reformulated in terms of nearest neighbours, so solution Eq. (4.4) of [TJ] applies. The 
appropriate substitutions are 

(m qL ,m^o,m lR ,m^o) -> {m- g ,m qL ,m^m^m^Jrn lR ) (2.4) 

As was discussed for the squark chain, primary distributions based on distinguishing 
the two leptons can in practice not be constructed. Instead of the distributions m qn i n , m qn i f , 
m qqln and m qql{ one can construct m qnl{high) , m qnl{low) , m qql{high) and m qql{low) . While the 
endpoint of a 'high' distribution is simply the maximum of the two primary endpoints for 
the given mass scenario, further considerations are needed for the 'low' distributions. For 
m q n lQow) these considerations are quite manageable, but for the three-particle distribution 
m qqi(\aw) they become rather involved. The general solution for the latter is given in App. [B]. 

The endpoint expressions are summarised below. Together with Eq. (|2.1|), the six first 
ones, Eqs. ( |2.5| )-( 2~T0| ), correspond to the distributions which are usable in that they do 



not rely on distinguishing the two leptons: 

KTf = H - O « - ™|o)/< ( 2 - 5 ) 

™^fiow) : see App. @ (2.6) 

maxKX,m^ f x ) (2.7) 



max 
rn qql{high) 
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(2.9) 
(2.10) 
(2.11) 
(2.12) 

(2.13) 
(2.14) 



m ™qk : E q- ( 4 - 4 ) of @ with ( m gL> m xO,m ifl ,m^o) -» (m§, m & , m^o, m^m^o/mjj 

(2.15) 



2.2 Quark ambiguity 

Adding a gluino to the cascade introduces an additional ambiguity: how does one distin- 
guish between the jet initiated by the gluino decay and that initiated by the squark decay? 
Initially we will assume that the two correct jets of a signal event have been selected, al- 
though we do not know which is which. (The difficulties associated with this assumption 
will be commented on later.) In this case, one can seek to distinguish the two jets kine- 
matically, by, for example, their transverse momenta. If m| — to| l is very different from 
m d ~ m -o> the two quarks will be emitted with noticeably different energies, and a pt 
comparison will, with a certain purity, serve to distinguish between the jets. However, one 
cannot, a priori, know which jet is to be assigned to q n or qt , i.e. whether it is q Q or which 
gives rise to the hardest jet, and one must instead turn to the differences in production 
rates. Only when a gluino sits on top of the chain is q n produced, and often, the ratio of 
directly versus indirectly produced squarks is sufficient to determine the average behaviour 
of qf . This then allows statements on the average behaviour of q n and q{ to be made, which 
may be used to construct samples where the quark identity is required, Eqs. ( |2.8] )-( ^To| ), 
as well as the squark distributions m qt u, m 9f ^ low ) and m qt i( high y Such a separation of q n 
and qf will always introduce impurities; whether these are manageable or not will depend 
on the given SUSY scenario. 

To be fully general one might consider constructing new secondary distributions in 
which there is no need to distinguish the jets. The distributions m;;, m qq , m qq u, m gg ^ low ) 
and mqq^high) would remain unaffected, but instead of m qf u and m qn u, one should con- 
struct new distributions m g ^(hi g h) and m q n^ ow y The endpoint expression of the former is 
readily available, but not of the latter. The complexity of the calculation and the resulting 
expression for w,^^ ow ) can be expected to be similar to that of w^fn ow ) 



shown in App. ||. 
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Furthermore, for the four m q i distributions one can collect four new distributions based 
on the order of the four values per event. Only the endpoint of the m 9 ^hi g h) distribution 
is readily available, while the other three endpoints become very difficult to calculate. A 
programme to calculate the new endpoints where quarks are not distinguished, has so far 
not been started. 

The above discussion shows that as a second jet is added, the complexity of the situa- 
tion increases significantly. In return many more endpoints become available, allowing for 
thorough consistency checks. In Sect. |3| we will return to these issues in the context of a 
simulation at the standard SPS la point. 

2.3 Theory distributions 

In order for the endpoint method to be useful, not only must the analytic expressions for 
the endpoints be available, it must also be possible to determine these endpoints from the 
experimental mass distributions. A first criterion for this is that the edges of the distri- 
butions 'point' unambiguously towards the exact endpoints. If the shape of a distribution 
is sufficiently concave at high values, the endpoint will most likely be underestimated and 
large systematic uncertainties must be added. In general the shapes of the mass distribu- 
tions vary if the sparticle masses are varied. To investigate the range of possible shapes 
for each distribution, the full cascade was generated for approximately 500 unique mass 
scenarios, and the resulting theoretical mass distributions studied visually. 

A representative selection of mass scenarios, showing some of the shape variety of the 
gluino distributions, is given in Fig. |l|. In generating these decays, quarks and leptons are 
assumed to be massless, which in the worst case (6-quark) leads to a perfectly acceptable 
error of O(MeV). Furthermore, only decay phase space has been used, without the associ- 
ated matrix elements. In principle, this ignores possible spin correlations between leptons 
and/or quarks. However, since the distributions used here are summed over lepton charge, 
no spin effects are expected [16]. A similar investigation for the squark distributions was 



performed in |12[j . See also |Tj 

Since m qq is constructed from two nearest neighbours in the decay chain, the shape is 
triangular for any masses, similar to run. Its endpoint can thus in principle be determined 
quite accurately. The shape of the m qq u distribution depends on the masses, but in prac- 
tically all scenarios the edge is well described by a linear descent towards the theoretical 
endpoint. This should guarantee its determination experimentally if the background is 
sufficiently low. For most scenarios this is also true for m qq i^ ow ^ , although the distribution 
can take on other shapes. 

The three 'high' distributions, m gg ;( higll ), m gf ^hi g h) and nt gjl z(high)> show a large variety 
of forms. The available shapes of the w- gf z(high) distribution are very similar to those of 
the ?n 9n z(high) distribution and are not reproduced here (see |12[| ). A common feature of 
these distributions is the danger of not noticing the 'foot', seen, for example, in m 9g z(high)> 
scenario (iv), where the true maximum value may be obscured by backgrounds. 

The last of the three-particle distributions, m qn u, also has an extensive variety of 
shapes. The behaviour of the edge is usually reasonable, but with a certain danger of foot- 
like structures, as in scenario (i). Finally, the two m„ B i distributions are often unreliable 
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Figure 1: Theory distributions for SPS la and four other mass scenarios, showing some of 
the shape variety. Only gluino distributions are shown. Kinematic endpoints are marked 
with a triangle and given in units of m™ff, the largest of the endpoints. (All pan- 
els start at zero invariant mass.) The scenarios (i)-(iv) are defined by the following sets 
of mass ratios r = (to^/to^ , /m^o, m^o/m^, m^/m^o): rp) = (1.03,1.30,1.07,6.01), 
r (ii) = (1.79,3.01,1.15,1.59), r (m) = (1.45, 1.05, 1.27, 5.08) and r (iv) = (1.03,1.29,2.02,3.60). 



in that the edges are concave. It may therefore be difficult to get good estimates of the 
endpoints. 



- 8 - 




Figure 2: SPS la cascade decay with branching ratios and cross-sections. 

In summary, from the 500 mass scenarios of which the distributions in Fig. |l| com- 
prise a fairly representative selection, the following can be stated regarding their overall 
behaviour: The distributions m qq and m qq u point reliably towards the endpoint for any 
mass scenario. Somewhat less reliable are the endpoint estimates found from the three- 
particle distributions. Lastly, even less reliable endpoint estimates can be obtained from 
the two m qn i distributions. In the next section we will see how these statements change for 
a 'realistic' experimental setup. 

3. Measuring gluino endpoints for SPS la 

Here we investigate the measurement of the endpoints associated with gluino decays for 
SPS la at ATLAS. Unless otherwise stated, the details of the analysis are identical to those 
for the squark chain, as described in Ref. |r^|, and will not be reproduced here. 



3.1 Signal and backgrounds 

The SPS la point is defined by the mSUGRA GUT-scale parameter values 

mi/2 = 250 GeV, m = -Aq = 100 GeV, tan/3 = 10, fi > (3.1) 



evolved down to the electroweak scale by version 7.58 of ISA JET [18]. The masses of particles 
at the electroweak scale can be found in fi2|[ . Fig. § shows the signal chain with branching 
ratios: to the left the production rate of gluinos is shown, followed by the branching ratios 
of a gluino into left-handed squarks and both sbottom states. Since the mass difference 
between b\ and the other squarks is comparable to the mass difference between the gluino 
and the squarks, phase space effects become significant, resulting in a noticeably enhanced 
decay rate into b\. 

Signal and background are generated for 300 fb _1 , which corresponds to three years at 
design luminosity of 10 34 cm _2 s~ 1 . The low-energy parameters are passed, via the standard 
interface, to PYTHIA 6.2 Jl9| which calculates the decay widths and the LHC cross-sections 



by use of CTEQ 5L |§ PDF's, and generates the events. Finally ATLFAST 2.60 |T| 
performs a parametrised fast simulation of the ATLAS detector. The simulation setup 
is identical to that used for the squark-endpoint analysis, except that four jets are now 
required with > 150, 100, 50, 20 GeV, rather than only three. Also, the definition of 
M e g is extended to include the fourth jet. 



- 9 - 



The QCD background is cut away by the requirement of two leptons 2 and of consid- 
erable missing pr- For the processes involving Z and W the requirement of high hadronic 
activity together with the missing removes nearly all events. After these rather hard 
cuts, the Standard Model background consists of approximately 95% ti. Most likely the 
Standard Model background is significantly underestimated, as the parton shower approach 
of PYTHIA does not generate at a realistic rate the hard jets which are required in the 
selection cuts. However, because of the large SUSY cross-section for this scenario, the main 
background will come from other SUSY processes. The possible underestimation of the 
Standard Model background is therefore not a danger for the current analysis. 

It is convenient to divide the background into three different parts, 'lepton-uncorrelated', 
'lepton-cor related' and 'combinatorial'. The two first consist of events which do not contain 
the signal chain. In the lepton-uncorrelated background two leptons are produced, but in 
different parts of the decay and therefore independently. With the assumption of lepton 
universality and neglecting the mass difference between electrons and muons, which is very 
reasonable for the energies involved, this background can be removed by subtracting the 
distribution for different-flavour leptons, as was done in Ref. |T^| . 

For the lepton- correlated background, the leptons are always same-flavour, so no 
different-flavour sample is available to show its distribution. These events typically come 
from the decay of Z, in which case they can to some degree be controlled, or from sleptonic 
decay of neutralinos. In the latter case no particular signature, e.g. in the mu distribu- 
tion, is available to discriminate this background from the signal. In particular much of this 
background will come from x^'s decaying sleptonically, but which are not part of our signal 
chain. For the analysis of the squark chain, as done in p^| , one important criterion is that 
the X2' 13 produced from squarks make up a good fraction of the total X2 production. Events 
with X2' 13 n °t originating from the relevant squarks usually still contain jets and constitute 
a large part of the background. As the signal chain grows longer, this background increases 
relative to the signal. In the case of the whole gluino chain what is important is that the 
x[>'s with a gluino grandparent are not severely outnumbered by the total \\ production, 
for the events selected. This is a stricter criterion, and we will find that the flavour of the 
intermediate squark, whether it is <jl or b, will be crucial for the isolation of the signal. 

The third background type, combinatorial background, is a result of our inability to 
know which of the jets correspond to the quark sent out from the gluino (q n ) and from 
the squark (<#). In most analyses one assumes that the quark from a squark decay can on 
average be distinguished from the quark from a gluino decay on the basis of hardness. In 
most mSUGRA scenarios this is viable, as X2 * s so much lighter than the squark. Since two 
squarks are expected in nearly all events, assuming here m g > mq, one can attribute the 
two hardest jets to the decay of the two squarks. Remaining jets can then be attributed 
to gluino decays, other quarks in the event, e.g. from the decay of t or W, and/or the 
underlying event. The exact jet selection algorithms used in this study are based on these 
assumptions, and will be detailed in the subsections below. 

2 For convenience of notation we will for the rest of the paper use 'lepton'/'slepton' for the two first 
generations and 'tau'/'stau' for the third generation. 
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Figure 3: Invariant mass distributions for non-6-tagged events. Endpoints not detectable. See the 
text for details. 



In more general SUSY scenarios it need not be the case that is harder than q n . In 
general, given an unknown SUSY scenario, the appropriate jet selection procedure must be 
the result of careful study. For this also the mass hierarchy of the gluino and the squarks 
must be established. A systematic study of how such information can be obtained is lacking 
but would be very valuable. 

3.2 Non-6-tagged distributions 

In Fig. | two of the non-6-tagged distributions, m qq and m qq u, are shown. In black with 
error bars the different-flavour-subtracted distribution ('SF-DF') is shown. Solid green 
shows the SUSY background ('SUSY'). Its shape is given by the lepton-correlated part, 
but the lepton-uncorrelated part is also significant and increases the random fluctuations. 
The Standard Model background is negligible and is not shown separately. The blue 
curve shows the parton-level distribution of the selected events ('TH'), while the red points 
with error bars show the part of the different-flavour-subtracted sample which contains 
the correct signal chain ('SC'), also referred to as the 'signal-chain distribution'. Any 
discrepancy between these two distributions is mostly due to combinatorial background, 
i.e. from picking the wrong jets. 

For these distributions was assumed to be one of the two hardest jets while q Q 
was selected among number three or four in pj--hardness. This is in line with the mass 
assumptions usually valid in mSUGRA scenarios, as described above. The jet selection 
used for the distributions plotted was, on an event-by-event basis, the one out of four 
possible combinations which gave the smallest m qq u value. In addition we required that 



neither of the involved jets was 6-tagged. (The 6-tagging simulation is described in [12] 
and is rather crude, although with a reasonable overall behaviour.) 

It is clear that the positions of the endpoints, which are seen in the parton-level distri- 
butions ('TH'), are not easily detectable from the different-flavour-subtracted sample. This 
is mainly due to the size of the lepton-correlated part of the SUSY background ('SUSY'), 



- 11 - 



which makes up ~ 80% of the total sample. Structures in the signal part of the sample 
are therefore not easily identified. This result can be anticipated from Table 3 of Ref. Jl2| , 
which shows that only 8.2/32.8 = 25% of q^s originate from a gluino. The ratio of gluino- 
induced x[>'s becomes therefore quite small. 

Another difficulty is the combinatorial background. The signal-chain distributions 
('SC') do not at all point to the nominal endpoints at 242 GeV and 490 GeV (for ul). 
This is because the jet pair selected is often not the correct one. 

Other jet selection algorithms of this simple type have been tried, but none allows any 
edge structure caused by the kinematics of the decay chain to be recognised. This is also 
true for the other five gluino distributions (not shown). One must therefore conclude that 
the endpoints of the non-6-tagged gluino distributions are not experimentally obtainable 
by looking at one distribution at a time. If instead correlations between mass distributions 
were investigated it might be possible to identify endpoint-related edge structures. 

3.3 6-tagged distributions 

The distributions of the 6-tagged samples are shown in Figs. |J-|5|. The curves follow the 
colour code of Fig. [3|, except the Standard Model background ('SM') is also shown (green 
points with error bars). For the 6-tagged distributions the different-flavour distribution 
is ~ 40% of the same-flavour distribution, which is somewhat larger than for the non-6- 
tagged distribution where the ratio is ~ 25%. In both cases the different-flavour distribution 
favours lower mass values and does not interfere very much with the edge structure. 

Contrary to the previous case, the 6-tagged distributions have clear edge structures 
which provide values for the endpoints. The main reason for this is that the different- 
flavour-subtracted SUSY background (solid green) now makes up a manageable 35%, to 
be compared with 80% for the non-6-tagged sample. This reduction is due to the fact 
that the majority of b 's are produced indirectly from gluino decay [because of the low b- 
content of the proton] in combination with a jet selection requirement of exactly two 
6-tagged jets. Although the total production of ^2 1S 6 - 7 times larger (see Fig. ^) than 
the production which starts out from g — > bb, two 6-jets are also needed, which reduces 
considerably the number of selected background events. (For background events the 6's 
are usually produced in the other chain, typically from the same g — > bb which then does 
not continue sleptonically, or from g — > t\t in either of the chains, which also produces 
multiple 6's.) Since the rate of 6-jets is considerably smaller than the rate of light jets, 
also in SUSY events, we are likely not to have additional 6-jets in a signal event. The 
combinatorial background is therefore small, as is clear from the good correspondence 
between the parton-level distribution ('TH') and the signal-chain distribution ('SC'). 

The specific jet-selection used here is in line with the previous assumption that 6f is 
harder than b n . The first is searched for among the two pr-hardest, the second among 
number three and four. Only events which have one 6-jet among the two hardest and one 
among the two next were selected. (In a more realistic study where emphasis is put on 
issues like fitting techniques, impact on the distributions from the precuts etc., it would be 
natural to investigate the effect of also including events which have their two 6-tagged jets 
as number 1 and 2 or as number 3 and 4.) 
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Figure 4: Four good 6-tagged invariant mass distributions. See the text for details. 



For the distributions in Fig. |5| no endpoint estimation is attempted. The two m^i 
distributions are distrusted since they do not usually 'point' towards the maximum value, 
as shown in Fig. |l|. The phase space corresponding to the largest values is very small, 
leaving only an experimentally undetectable tail to mark the endpoint. These distributions 
are furthermore made less useful by the need to distinguish b n from bf . For the distributions 
of Fig. | this is not an issue because both 6's were used. At SPS la the endpoints involving 
6 n have smaller nominal value than the corresponding endpoints with bf. Consequently, 
if we are not able to cleanly distinguish b n from 6f, all distributions which have b Q as 
the only quark will get a contamination which stretches beyond the endpoints. If this 
contamination is substantial, and here it is, the b n endpoints will be washed out. The 
combination of these two effects seems to disable these three distributions. The nominal 
endpoints are m^f = 281 GeV, m™f low) = 197 GeV and m^f Mgh) = 248 GeV, which 
are clearly not obtainable from the distributions in Fig. ||. To cope with the problem of 
distinguishing b n and 6f more general distributions could be constructed, as was discussed 
in Sect. p.2[ 
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Figure 5: Three bad 6-tagged invariant mass distributions. See the text for details. 



For mfcb the fit values are a considerable 25 GeV higher than the nominal value, and 
are not really accounted for by the errors. It is not completely understood why the dis- 
crepancy is so large; whether it is a statistical fluctuation or some other effect. Since the 
reconstruction resolution is much worse for jets than for leptons, one should expect that 
a mass distribution constructed from jets alone will be more smeared and therefore give 
worse endpoint determination than one which also involves leptons. This is probably part 
of the reason. 

To extract the endpoints from the distribution, signal and background hypotheses are 
needed. Usually, at this level of detail, a straight line is used for the signal in the edge 
region. We see from the theory distribution of Fig. [I] that the straight-line hypothesis in 
the edge region is normally quite well supported for the four distributions of Fig. |4|. (The 
exception is rn qq i(high) i which in some cases might have a dangerous foot at the very last 
part of the edge.) 

Which background hypothesis to use is less clear. The background is here mainly from 
other SUSY processes, and therefore in principle unknown. One attempt to sketch the 
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shape of the background, is by combining the lepton and the jet sectors of different events 
into a mixed-event sample, as described in Sect. 5 of [p^ |. The shapes of the mixed-event 
distributions are then used as background hypotheses. Another way is to simply select for 
the background some appropriate function, based on the distribution slightly beyond the 
edge, where only background resides. An exponential or a polynomial usually gives a fair 
description. Whatever the signal and background hypotheses, systematics are introduced 
which may not always be easy to estimate. 

For the plots of Fig. ^| the signal was fitted with a straight line, while the background 
was modeled by both an exponential and a mixed-event sample. A four /three-parameter 



fit was then performed by use of MINUIT [22] for exponential/mixed-event background 



hypothesis and several different histogram binnings and fit ranges. Each distribution gives 
fit values and statistical errors within a fairly narrow interval, indicative of the systematic 
uncertainty in the fitting procedure. The results of the fits are summarised in Table |l] at 
the end of this section, after the discussion of energy scale errors. 

3.4 Propagation of energy scale errors 

In ATLAS it is expected that the absolute energy scale of jets and leptons will be known to 



1% and 0.1%, respectively, see Ch. 12 of [23|. This energy scale uncertainty translates into 
an uncertainty, or 'error', for any mass constructed from jets and/or leptons. In particular, 
the masses which go into our distributions have such an error in addition to the statistical 
error. 

Invariant masses which are constructed from jets alone or leptons alone, will inherit 
their uncertainties of 1% or 0.1%, respectively. For an invariant mass made from one lepton 
and one jet, the energy scale error is at 0.5%. For masses constructed from more than two 
jets and leptons, the energy scale error is not constant, but varies within a given calculable 
interval, depending on whether it is the jet(s) or the lepton (s) which dominate the invariant 
mass for the given event. 

To find appropriate energy scale errors, we must investigate m qq u, m qq i and m q n. The 
overall masses can be expressed in terms of two-particle masses, which have a constant 
energy scale error: 

m \ q ii = m \q + m2 gn i n + m2 q j { + m \i^ + m l ( i ( + m a ( 3 - 2 ) 

m lql = m lq + m L + m ll ( 3 - 3 ) 
m \ll = m2 ql n + m2 ql ( + m l ( 3 - 4 ) 

Absolute lower /upper limits for the energy scale error of the quantities on the left-hand 
side can then be found by assuming that the right-hand side is totally dominated by the 
term which has the smallest /largest energy scale error. This results in 

ft^e (o.i,i)%, (0.5,1)%, (0.1,0.5)% (3.5) 

m qq u m qq i m q ii 

where a denotes the energy scale error. These are absolute limits valid for any mass 
scenario. To find the relevant numbers for SPS la, all accepted events were reexamined; 
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scaling the jet momenta of each event by 1.01 and taking the ratio of the new and the 
old invariant mass. We then find the following average and root-mean-square values (in 
parentheses) of the relevant energy scale errors, 



^("^(high)) 
m M>i(high) 



0.66(0.10)%, 
0.71(0.10)%, 



f (»™NW(1ow)) 
(low) 

°( m b n ll) 



0.80(0.11)% 
0.42(0.02)% 



(3.6) 



Inclusion of the lepton energy scale will give a small correction to these numbers. For 
m b u ih which we will not be using, the energy scale is nearly constant. For the three other 
distributions, the errors lie between approximately 0.5% and 1%, and fairly uniformly 
distributed, as is reflected in the root-mean-square values. For the fitting, only the energy 
scale error for masses which lie in the edge region is relevant. However, it turns out that 
at SPS la and for these distributions, the error is almost the same for low as for high 
invariant mass values. Although the error of each distribution is found to lie in a fairly 
broad interval rather than being constant for all events, we have here used the average 
values, Eq. ( |3.6|) , as a basis for the energy scale errors in Table @. 

An alternative approach could be to scale the jet momenta up/down as done above, 
then redo the entire fitting process on the new distributions and from this extract the effect 
of the energy scale error. For this to work one would have to disentangle the effect of the 
scaling from the yet uncontrolled systematics of the fitting procedure. At the present stage 
of fitting competence the gain from using a more correct procedure is probably lost in the 
increased complexity. 





Nominal 


Fit 


Energy Scale 


Statistical 


Edge 


Value 


Value 


Error (cr scale ) 


Error (cr stat ) 




[GeV] 


[GeV] 


[GeV] 


[GeV] 


m bb 


312.7 


335-339 


3.4 


6-10 


m£Lx 
m bbll 


496.3 


494-500 


3.3 


5-7 


iyy} ITleLX 

''*&&Z(low) 


413.2 


407-417 


3.3 


8-12 


/YY} fflelX 

" 4 M>2(high) 


461.9 


454-462 


3.3 


5-7 



Table 1: Endpoint values found from fitting the edges in Fig. [|, for 300 fb 1 . The nominal values 
correspond to the mass of b\ which is produced at significantly higher rates than the heavier &2- 



The energy scale errors are based on the discussion in Sect. 3.4. No values are given for the three 
distributions in Fig. ^| 



4. Masses from endpoints 
4.1 10,000 ATLAS experiments 

The sparticle masses can be obtained from a numerical fit based on the gluino endpoints 
found in the previous section together with the squark endpoints of [O]. Due to the 
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simplistic treatment of the detector effects and the endpoint fitting, we do not trust the 
fit values to be representative of what a more complete simulation would give. Therefore, 
rather than use directly the endpoint values obtained from the distributions, we generate 
an ensemble of 10,000 ATLAS experiments based on the estimated errors, as done in |10| . 

For each experiment the endpoint values E exp are combined with the general endpoint 
expressions E th in a least-square function E, 

E = [E cxp - E th (m)] T W[E cxp - E th (m)] (4.1) 

The weight matrix W is the inverse of the covariance matrix which is constructed from 
the endpoint errors and appropriately handles the endpoint correlations due to the energy 
scale error. The minimisation of E then yields the masses. Due to the composite nature 
of the endpoint expressions, there are usually several competing E minima for a given set 
of endpoints. If these minima are close in S value, they must all be considered, opening 
the way for multiple mass solutions. Finally, ensemble distributions for the masses can be 
plotted and studied, and are to be interpreted as probability distributions for the masses 
which can be obtained at the LHC. We will be interested in the mean and width of these 



ensemble distributions. For more details of the procedure, see [12|. 
4.2 Mass estimation 

In total six sparticles are involved in our decay chains: Xi> Iri X2> QL, bi and g. The 
masses of the first five of these could already be obtained from the squark endpoints alone, 
as was done in |l2]. The gluino mass is accessed by the gluino endpoints. Since these also 
involve the other masses (except rriq L ), the squark endpoints and the gluino endpoints are 
combined in a single fit where all the masses are free parameters. The six squark-endpoint 
measurements are given in Table 4 (upper part) of The values of the right-most 

column, 'Syst. Fit Error' are not used. For the gluino endpoints of Table [l] (in the present 
paper) the statistical errors are given as intervals. We take here the midpoints of the 
intervals, and further assume that the systematics of the fit values (also given by intervals) 
will be settled, or at least dominated by the statistical errors. This assumption might be 
somewhat optimistic, and should be kept in mind when contemplating the results obtained 
below. (In the next section we will look at the effects of doubling the errors.) Furthermore, 
we decided to use only three out of seven gluino distributions. In particular the choice of 
excluding ?n^ ax was made after consultation with the Monte Carlo truth at our specific 
SUSY scenario, which is not a viable strategy in a realistic setting. More study could 
however promote this result into a generic one. 

If we start from the squark endpoints and add only one of the gluino endpoints in 
the numerical fit, the gluino mass will be determined and all the other masses will remain 
unchanged. This is so because there is exactly one new measurement for one new mass. 
The gluino mass returned from the fit is the one which gives zero for the new terms added to 
E, so there is no increase in the E-value. Also, the number of minima remains unchanged 
by adding one more measurement for one more mass. This situation is similar to the 
determination of mg i in [12]. Only one endpoint involves &i, so its inclusion has no effect 
on the other masses or number of minima. 
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When more gluino endpoints are added, the gluino sector becomes over-determined, 
and the positions of the minima will change, i.e. the other masses will be affected. Since 
the gluino endpoints have somewhat larger errors than the squark endpoints, large effects 
are not expected, except perhaps for b\. Below, results are given for the case when all three 
selected gluino endpoints are used. 

For the current precision of the endpoint 
measurements, the numerical fit nearly always 
returns two minima, one in mass region (1,1), 
which is the region of the nominal masses at 
SPS la, and one in mass region (1,2). If the 
minima are close in S-value, both must be 
considered. Table || shows the probability of 
having more than one solution in an exper- 
iment. The cut on AS, the distance to the 

global minimum, gives the quality of the second minimum. In most cases the (1,1) mini- 
mum is the selected one. These numbers are very similar to the numbers obtained without 





# Minima 


(1,1) 


(1,2) 


AS = 


1.00 


91% 


9% 


AS < 1 


1.11 


95% 


16% 


AS < 3 


1.30 


98% 


32% 


AS < 99 


1.87 


99% 


87% 



Table 2: Number of minima for various AS 
cuts and their whereabouts. 



the gluino endpoints, Table 5 of [12], where a more detailed description can also be found. 







(1,1) 


(1,2) 




Nom 


Mean 


RMS 


Mean 


RMS 


Al 


96.1 


96.3 


3.7 


85.5 


3.4 


m ln 


143.0 


143.2 


3.7 


130.6 


3.8 


A2 


176.8 


177.0 


3.6 


165.7 


3.5 




537.2 


537.5 


6.0 


523.5 


5.0 


m ' bl 


491.9 


492.2 


12.5 


471.8 


12.6 


m-g 


595.2 


595.5 


7.2 


582.5 


6.8 


lR Xi 


46.92 


46.93 


0.28 


45.11 


0.72 


m^o — m^o 

A 2 Xl 


80.77 


80.77 


0.18 


80.19 


0.29 


mq L - m^o 


441.2 


441.2 


3.1 


438.0 


2.8 


m hi - m^o 


395.9 


396.0 


11.2 


386.3 


11.2 


m g - m^o 


499.1 


499.2 


5.6 


497.0 


5.4 


m-g - m 5i 


103.3 


103.3 


9.1 


110.7 


9.5 



Table 3: Masses (Mean) and root-mean-square deviations from the mean (RMS) of minima in 
regions (1,1) and (1,2), for AE < 1. The nominal masses (Nom) are given in the second column. 
All values in GeV. See the text for more details. 



The masses are given in Table || for minima which satisfy AS < 1. The numbers are 
very close to the ones obtained without the gluino endpoints, see Table 6 of |l2|]. Only b\ 
is affected, as was expected. 

The gluino mass is quite well determined. The ensemble mean is at the nominal value, 
and the root-mean-square deviation from the mean is 7.2 GeV, only a GeV more than for 
For the other masses, especially the lighter ones, mass differences are much more 
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accurately determined than the ma sses themselves. This is also the case for the gluino, 
3 The definition of the regions is given in Sect. 4.3 of For a mass region the first index refers to 
the expression used for mJJJff, the second index shows which combination is used for (jri^^(i ow ),m-^t(high))i 
see Eqs. (4.4)-(4.5) of Q . 
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Figure 6: Sparticle masses and mass differences for solutions with AS < 1. See the text for details. 



although to a lesser degree, as is seen from the smaller root-mean-square value for m g — m^o 
of 5.6/5.4 GeV. When the correlation to Xi is taken into account, as in mass differences, 
ijL is better determined than g, as it naively should be from the smaller squark endpoint 
errors. Since the gluino endpoints involve the sbottom mass, one might expect that the 
masses of g and b\ have a considerable correlation. Some correlation is found, as can be 
seen from comparing the root-mean-square value of m g — mg with that of alone. 

The ensemble distributions summarised in Table ||, i.e. for AS < 1, are shown in 
Fig. ||. The unfilled distributions in black show from left to right m^o, rrif , Tt^o, m^, 
rriq L and m g for solutions in the nominal region (1,1). We will have such a solution in 95% 
of the experiments, see Table |2[ The unfilled distributions in blue (grey) show the same 
masses for solutions in region ( 1,2). Such a solution occurs in 16% of the experiments. (For 
AS < 1 there is an 11% chance of having both solution types.) The smaller rate of the 
(1,2) solutions is reflected in the smaller area under the blue curves. The (1,2) solutions 
return slightly lower masses, so these distributions are positioned 'down' (due to the lower 
rate) and slightly to the left of the (1,1) distributions. The filled distributions show from 
left to right m^—m^o, m^o—m^o, rrig—m^, m^—m^o, rriq L —m^o and m§— m^o. Again, 
the most populated distributions (black curves) are for solutions in region (1,1), the least 
populated (blue curves) for (1,2) solutions. For mass differences there is more overlap 
between the (1,1) and (1,2) solutions, in particular for —m^o and m^o—m^o, of which 
only the lower parts of the distributions are visible. Mass differences are better determined 
than the masses themselves, reflected here by the narrower distributions of the former. In 
the case of b± this effect is clear for m g — , much less so for — m^o , which points to 
the fact that b\ largely decouples from x? ; 
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4.3 Edge sensitivities 



It may seem a puzzle that nig is determined with much higher precision than m^. Both 
masses enter the three new gluino endpoints, but has in addition the m ^ 1 if^g > ]L) mea- 
surement. Naively one might therefore expect the b\ mass to be more precisely determined. 

In order to better see how masses are constrained by the endpoint formulas, consider 
the part of the least-square function £ which involves one given endpoint and one given 
mass only. For simplicity of argument, assume that the endpoint measurement is at the 
nominal value, as are the masses. Changing the given mass away from the nominal value 
by an amount Am will give an increase of our selected S part by 

/E(m + Am) - £ ex P\2 ^ /E[m) + (dE/dm)Am - £ CX P\2 _ /(^/3m)Am\2 

where E is the value of a particular kinematical endpoint. The value of Am which gives 
an increase of £ by 1 (from this term only), is then given by 

Am = - (4.3) 

(dE/dm) v ; 
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0.6 
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m ?f«(high) 


-3.2 
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1.0 


-0.3 


-8.2 


0.8 


2.7 










min 
m b f ll(d>%) 


-0.7 


-6.5 


-1.5 


-3.1 


1.5 


3.0 






0.4 


11.8 






m bb 










-0.3 


-39.6 






-1.3 


-8.1 


1.7 


6.3 


" l bbll 


-1.3 


-5.2 






0.3 


23.4 






-0.1 


-65.1 


1.0 


6.5 


" l bbl(law) 






-2.3 


-4.5 


1.5 


6.8 






-0.6 


-19.0 


1.3 


8.4 


max 
'"■66i(high) 


-2.1 


-3.2 


1.4 


4.7 


-0.4 


-17.6 






-0.3 


-24.5 


1.1 


6.1 


max 
m b n ll 

yy^fflcLX 

6 n i(low) 

™ IIlclX 

m 6„/(high) 


-1.2 
-2.1 


-2.6 
1.4 


0.7 
2.1 




-1.2 
-0.9 
-1.1 


1.5 
1.0 
1.3 



Table 4: Partial derivatives of endpoints with respect to masses at the nominal mass values. The 
larger the partial derivative dE/dm is, the more sensitive the given endpoint is to the given mass. 
The quantity Am [GeV] is defined by <te / '(dE / 'dm), where ge is the combined statistical and 
energy scale error of the endpoint. Numbers are specific for SPS la. See the text for more details. 

Since the shift in mass also induces changes to other contributions to S, the inter- 
pretation of Am is not straight-forward in terms of mass error. Still, Am does combine 
the precision of the endpoint, given by the endpoint error ge, with the response of the 
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endpoint to the mass, which are the two important quantities. The same relation could be 
derived in a mathematically more intuitive way, 

where the endpoint variation AE is set equal to the experimental error erg. 

Inclusion of the endpoint error opens for comparison between different endpoints, even 
if only in a semi-quantitative way which does not take correlations into account. The 
importance of the different endpoints in constraining a given mass can then be studied. In 
Table [|, values for dE/dm and Am are given for both the squark and the gluino endpoints. 
For the three gluino endpoints of Fig ||[ where no fit was made and hence no a is available, 
Am is not given. 

These numbers show immediately that the three gluino endpoints included in the 
numerical fit have a strong sensitivity to the gluino mass, but only a very moderate sen- 
sitivity to the sbottom mass. Also wij^^s) is seen not to constrain too much. If 
m^ ax could be used or some of the squark endpoints other than ^j^^e) (numbers shown 

only for t^^"^-))' then might be determined more accurately. The table can also 
be used to understand better which of the endpoints are important in constraining each of 
the lighter masses. 



5. Increasing the errors 

In the previous section the uncertainties of sparticle masses were found from the esti- 
mated errors of the endpoint measurements. The technique of generating an ensemble of 
experiments with endpoint measurements spread Gaussian-like around the nominal values 
sidesteps some of the problems of the crude detector simulation and edge-fitting procedure, 
as only the estimated statistical errors and energy scale errors are then relevant. However, 
these are also subject to uncertainties. First, for the statistical errors we used the means 
of supposedly reasonable intervals obtained from recurrent fits with slight variations in bin 
width, range of fit and fit functions. Clearly these numbers are tentative. Second, one 
cannot be certain that the energy scale will actually be known to the precision aimed for. 

In this section we will therefore investigate two somewhat more conservative scenarios. 
First we double the statistical errors of all the endpoints, keeping the energy scale errors at 
the 'default', i.e. at the values given in Table [j] (this paper) and Table 4 of |fl3| . Then we 
keep the statistical errors at the default values and double instead the energy scale errors. 

However, in order to better understand the changes this will lead to, we first try to 
disentangle the impact of the statistical errors from that of the energy scale errors for a 
simplified version of the default scenario considered in the previous section. 

5.1 Disentangling statistical from energy scale errors 

A simplified method is considered, in which only the four endpoints m™ ax , m™f*, m ^^ ow ) 

and "z™f(high) are used to obtain the masses of Xii Iri X2 ari d Ql- The main features 
of the comprehensive method where all the endpoints are used, are also present for this 
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simplified analysis. In particular we can use the inversion formulas of [p^| . For region 
(1,1) Eqs. (4.13)-(4.16) are appropriate. Notice in these formulas how the three lightest 
sparticle masses are proportional to m™ ax . Hence, these three masses scale with the lepton 
energy scale error. Conversely, compare this with their dependence on the three squark 
endpoints. If these three endpoints have identical relative energy scale errors [which in 
reality is only approximately true for m™f* - here we have assumed it to be exactly true], 
then the energy scale errors cancel out for the three lightest masses, as is easily seen by 
letting each of the squark endpoints be scaled by a factor (1 + 5). Consequently, the impact 
of the energy scale errors on the three lightest sparticles is equal to the lepton energy scale 
error, which for most purposes is negligible. For the squark mass, this symmetry is broken, 
rriq L is not proportional to m™ ax , so the other relations do not apply either. The squark 
mass therefore has a noticeable dependence on the energy scale errors. 

In order to quantify the statements just made, 
an ensemble of ATLAS experiments was con- 
structed in the usual way [12], using the de- 



fault errors of the four squark endpoints. Ta- 
ble H shows the root-mean-square values thus 
obtained for the four masses and the mass dif- 
ferences. Three cases are shown: only statistical 
errors, only energy scale errors and combined. 
Even though the energy scale errors are roughly 
twice the statistical errors for the endpoints be- 





Only 


Only 


Both 




^.stat 


^.scale 




m x? 


3.9 


0.10 


3.9 


m ln 


3.9 


0.14 


3.9 


m x% 


3.8 


0.18 


3.8 


mq L 


5.9 


2.4 


6.3 


lR Xi 


0.28 


0.05 


0.28 


m^o — m^o 

X.2 Xl 


0.16 


0.08 


0.18 


rriq L — m^o 


2.0 


2.4 


3.1 



Table 5: Disentangling statistical from en- 
ergy scale errors. Root-mean-square values 
of uncertainties in masses and mass differ- 
ences [GcV]. Only the four basic endpoints 
have been used. 



ing used (Table 4 of [12]), the latter clearly dom- 
inate the mass errors. Only for combinations 
involving do the energy scale errors play a 
role. For mq L —m^,o the impact of the two uncer- 
tainties are comparable. Comparison with Ta- 
ble H shows that the stripped scenario accounts 

for the main features of the complete scenario to which we return below. 

5.2 Doubling the statistical errors 

Having established that the statistical er- 
rors dominate the errors on the masses, we ex- 
pect large effects from doubling these errors. 
Table |6| shows a noticeable increase in the rate 
of (1,2) solutions relative to the nominal (1,1) 
solutions. This is no surprise, as SPS la lies 
very close to the border. Allowing larger er- 
rors for the endpoints naturally allows more of 
the (1,2) region to be accessed. Furthermore, 
there is an increase in the number of solutions, 

which is also reasonable. The 10-15 GeV lower mass values of the (1,2) solutions become 
more troubling as the occurrence of this minimum becomes more frequent. 





# Minima 


(1,1) 


(1,2) 


AS = 


1.00 


73% 


27% 


AS < 1 


1.24 


84% 


41% 


AS < 3 


1.47 


90% 


57% 


AS < 99 


1.87 


93% 


75% 



Table 6: Number of minima for various AS 
cuts and their whereabouts. The statistical er- 
rors are doubled, compared to Table |[ 
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m l R 

A2 

nig L 



m iR 



in. 



<1L 



111 



bi 



- m^o 

Al 

- m^o 

- m^o 

Al 
Al 



Nom 



96.1 
143.0 
176.8 
537.2 
491.9 
595.2 



46.92 
80.77 
441.2 
395.9 
499.1 
103.3 



(1,1) 
Mean RMS 



97.5 
144.5 
178.3 
539.3 
494.4 
597.4 



7.2 
7.0 
7.0 
11.0 
24.3 
12.9 



47.03 
80.80 
441.8 
396.9 
499.9 
103.0 



0.53 
0.32 
4.4 
21.9 
9.4 
18.0 



(1,2) 
Mean RMS 



84.2 
129.0 
164.3 
522.1 
468.6 
581.1 



44.72 
80.07 
437.8 
384.4 
496.9 
112.5 



6.1 
6.7 
6.2 
8.7 
24.0 
11.7 



1.21 
0.48 
3.7 
21.9 
9.0 
18.9 



Table 7: Same as Table except for doubled statistical errors: Mean and root-mean-square 
deviations from the mean (RMS) [GeV]. 



From Table [?] we see that doubling the statistical errors very nearly amounts to doubling 
the errors on the masses. The exception is mainly rriq L — m^o whose error increases less. 
This effect is explained by the discussion in Sect. |5.1| , see Table [|, which shows that for 
this particular mass combination the statistical errors only account for half the total error 
on the mass difference. For g a pattern similar to, but less pronounced than that of qi can 
be seen. 

5.3 Doubling the energy scale errors 

As already anticipated by the previous discussions, changes in the energy scale errors have 
only minor effects in the SPS la scenario at the given integrated luminosity. The doubling 
of the energy scale errors does not at all change the number and positions of minima from 
the default values of Table ^. 

Furthermore, as is seen from Table ||, most masses remain largely unaffected, the 
exception being combinations involving m^ L and to a lesser extent rrig, as previously noted. 
The explanation is the same: For the lighter masses the jet energy scale errors cancel out, 
leaving only the minimal effect of the lepton scale error. In the case of niq L it does play a 
role, although not dominant, giving a slight increase, and for vaq L — m^o the impact of the 
energy scale errors is again significant. The same applies to m g and rrig — rn^o. 

6. LHC + Linear Collider (LC) 

Within the time-frame of the analysis of LHC data, measurements from a Linear Collider 
may become available. While the LHC is able to measure mass differences at high precision, 
as documented in the previous sections, a Linear Collider will, due to the much cleaner 
environment, provide precise mass measurements of the sparticles which are kinematically 
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(1,1) 


(1,2) 








RMS 




RMS 


TYl ~0 
Xi 

rtiT 

l R 


96.1 


96.3 


3.7 


85.5 


3.4 


143.0 


143.2 


3.7 


130.6 


3.8 


x 2 


176.8 


1 77 n 


3 6 


165 7 


3 5 


,n qL 




537.3 


7.3 


523.8 


6.4 


m- h 


491.9 


492.1 


13.0 


472.1 


13.1 


9 


595.2 


595.3 


9.3 


582.9 


8.9 


lR Xi 


46.92 


46.93 


0.29 


45.11 


0.72 


m x% ~ m x? 


80.77 


80.77 


0.22 


80.19 


0.32 




441.2 


441.1 


5.1 


438.3 


4.9 


m- bi - m^o 


395.9 


395.8 


11.8 


386.6 


11.8 


m- g - m^o 


499.1 


499.0 


8.1 


497.4 


8.1 


m g - m h 


103.3 


103.2 


9.3 


110.9 


9.8 



Table 8: Same as Table g], except for doubled energy scale errors: Mean and root-mean-square 
deviations from the mean (RMS) [GeV]. 



accessible, in particular Xi- Such a measurement will be the scale fixer which is lacking in 
the LHC data, and will in combination with the LHC measurements allow one to also fix 
the masses of the heavier sparticles not accessible at the Linear Collider. 

To estimate the effect of a Linear Collider measurement of m^o, we add to our least- 
square function £, a term [(m^o — ™~?)/<7m C ] 2 ' wnere the quantity with superscript 'LC is 
the Linear Collider measurement, and generate an ensemble of LHC and LC experiments. 
Since <7^ c = 0.05 GeV [24], the LC measurement practically fixes m^o at the nominal 

value. (For the LHC measurements we use in this section the default endpoint errors.) 

For the (1,2) solutions, which normally re- 
turn Xi masses some 10 GeV below the nom- 
inal value, the LC measurement has the dra- 
matic effect of reducing their occurrences to ~ 
1% (for AS < 3). These minima can there- 
fore for most purposes be neglected. As a con- 
sequence, the probability of having two min- 
ima is strongly reduced, to the per mille level. 
For SPS la the Linear Collider measurement 
thus closes the issue of multiple minima alto- 
gether. (This behaviour was already reported 
in fi"^| , the inclusion of the gluino endpoints 
thus makes no difference in this respect.) 

The combined LHC + LC results are shown 
in Table ||. Comparison with the numbers of Table || shows that the precision of the mass 
determination improves considerably when the LC measurement is included. The root- 





Nom 


(1,1) 
Mean RMS 


Xi 


96.05 


96.05 


0.05 


lllr 


142.97 


142.97 


0.29 


m^o 

X2 


176.82 


176.82 


0.17 


mg L 


537.2 


537.2 


2.5 


m ~ bl 


491.9 


491.9 


10.9 


m-g 


595.2 


595.2 


5.5 


m-g - m gi 


103.3 


103.3 


9.0 



Table 9: Mass values (all in GeV) from 
LHC+LC. Since the occurrences of (1,2) solu- 
tions are reduced to ~ 1%, they are left out. 
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mean-square values of the masses are now approximately equal to the root-mean-square 
values of mass differences without the LC measurement. For mr the correlation to m^o 

01 Xi 



m?, ~ mr 



is not dominant. The spread of is therefore still larger than the spread of 
which does not feel the fixing of m^o. As was the case without the LC measurement, the 

° Xi 7 

inclusion of the gluino endpoints in the numerical fit does not affect the other masses. The 
exception is m^, for which the root-mean-square value is reduced by ~1 GeV compared 



to the results in [12]. 
7. Conclusion 

In this paper we have extended the endpoint method of obtaining masses in R-parity 
conserving SUSY scenarios to also include the gluino mass, given the decay chain (|1.2|). 
All the gluino endpoints have been calculated. Details of the calculations of and 
m q l { tr are given i n the appendices. Theory distributions for the seven new distributions 
were studied for a selection of mass scenarios, revealing that some of the distributions often 
have little phase space towards higher masses, making them less useful. 

An ATLAS simulation of 300 fb _1 was performed for the mSUGRA point SPS la. 
While we were not able to detect the gluino edges in the case of an intermediate first or 
second-generation squark ((Jl), edges were found and fitted in the case of an intermediate 
sbottom. This is due to fact that 80% of the b's come from a gluino, together with the jet 
selection requirement of having two and only two 6-tagged jets, both important background- 
reducing factors. Not all of the endpoints were used. Energy scale errors were discussed 
and found to lie in intervals. Yet, they were taken as constants in the analysis. 

To estimate the precision with which sparticle masses can be obtained at the LHC, an 
ensemble of 10,000 'gedanken experiments' were performed, in which for each experiment 



three of the gluino endpoints were combined with the squark endpoints obtained in [12] 
in a least-square fit to give the masses. Mass differences are better determined than the 
masses themselves, a generic feature of the method. Furthermore, one set of endpoint 
measurements in general corresponds to several sets of masses. The inclusion of the gluino 
endpoints affects the number of minima and the masses of Xi> Ir, X2 an d Ql only minimally. 
In the case of b\ there is a noticeable correlation to the gluino. The ensemble distribution 
of the gluino mass was found to have a root-mean-square value of 7 GeV. The spread for 
rrig — m^o was found to be about 1.5 GeV smaller. 

Spurred by the lack of significant improvement in the ensemble spread of compared 
to the precision obtained for nig, even though the new measurements involved both, a 
sensitivity analysis for all the relevant endpoints was performed. This investigation showed 
that the three gluino endpoints being used have limited sensitivity to the sbottom mass, 
thus confirming the different impact they have on the gluino and sbottom masses. 

Two alternative error scenarios were investigated, one in which the statistical errors 
obtained from the fast simulation of SPS la were doubled, another in which the energy scale 
errors were doubled. It was shown that the statistical errors of the endpoints, although 
smaller in magnitude than the energy scale errors, dominate the error on the masses. Only 
in a few exceptions, like niq L — m^o, do the energy scale errors play a significant role. 



- 25 - 



Finally the impact of a joint LHC-LC analysis was estimated. The endpoint measure- 
ments from the LHC were combined with a Linear Collider measurement for the LSP mass. 
This essentially fixes the SUSY mass scale. Consequently, in the combined analysis the 
masses themselves are determined with roughly the same precision as that of mass differ- 
ences determined at the LHC alone. According to a dedicated study [25], the anticipated 



accuracy will suffice to determine the high-scale mass parameters with a precision ranging 
from 0.1% for m\j2 to 14% for Aq. 

A. Calculation of m^f 

Below we use the calculation of the known squark endpoint m™^ x as an aid to obtain m™^f. 
We first follow what may seem to be the most obvious route, then note that while this works 
fine for ?n™^ x , the increased complexity of having one more particle to consider strongly 
limits the applicability to m^f. An alternative method is then developed by which m^f 
can be found in a straightforward manner. 

A.l Method 1: By means of angles 

In the first method we keep the calculation close to the physical progression of the decay, 
picking up all decay angles and eventually maximising with respect to these. Fig. |7| shows 
in a stepwise manner the qi decay chain from which m™^ x is calculated. Start in the rest 
frame of qx, and align the coordinate system to its decay products. Next, boost to the rest 
frame of x®, which decays such that l n is emitted at an angle a relative to the direction of 
Qf. Finally, boost to the rest frame of Ir, which in turn decays; for a given angle a, m qt u 
is maximised by choosing p if opposite to p„ f + p in in the rest frame of Ir [giving a planar 
decay configuration]. 

The resulting m qi u is expressed in terms of the angle a and the four sparticle masses 
involved. In order to maximise with respect to a, the critical points of this constrained 
m qi ll distribution must be sought (and tested if they are maxima or minima). In the part 
of mass space where no critical point is found, the endpoints of the domain, a £ [0, tt] , must 
be analysed. In this case, three different expressions are found for three distinct regions 
in mass space. In total, mass space is thus divided into four regions, each with its specific 
expression for ?n™^ x , see Eq. (4.4) of [12]. (See [|ll]] for a different derivation.) 



The above calculation becomes fairly involved, and is preferably performed with the 
assistance of a computer program. If we now want to calculate rn^f, a gluino is added 
at the head of the chain. In general two more angles are needed to describe the gluino 
decay. However, since the maximum of the endpoint values are obtained with planar 
decay configurations, one of the angles can be dropped. The remaining additional angle is 
nevertheless sufficient to complicate significantly the search for extrema of m qq u . 

A. 2 Method 2: On a line 

In the alternative approach to m™f x (and m™^) the physical progression of the decay is not 
tracked. Instead an expression for the endpoint is found in terms of the one missing four- 
vector, that of X?- Arguments which invoke the notion of dominant decays are then applied 
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Figure 7: Traditional way of calculating m™ff . 

to point out a number of limiting decay configurations (of great geometrical simplicity) 
which must be considered. 

Investigating first m^ x , the four- vector of the sum of the end products; the two 
leptons, the quark and the LSP, equals the four-vector of the squark. We therefore have 

m \ L = (PqtU + Px?) 2 = m \ t U + m \\ + 2 Pq t ll ■ Pxl ( AA ) 

where p qf u is the sum of the four-vectors of the two leptons and the (far) quark. In the 
rest frame of Ijl we have 

PqtU = ~Px?' E 9tU = m 9L ~ E x1 ( A - 2 ) 

which gives 

m lll = m L ~ m \\ ~ 2 { m qL\] m \o+V>\a - m |o) (A.3) 

The smaller the momentum of the LSP is in the rest frame of the initial squark, the larger 
m„ { n is, which agrees with one's intuition. At the extreme, if the LSP is brought to rest in 
the squark rest frame, m q{ u will attain its largest value, rriq L — rn^o, which corresponds to 
the critical-point solution of method 1. 

In some regions of mass space it is however not possible to have Xi a t rest in the rest 
frame of qi- In these regions, one of the intermediate sparticles is sent off with so high 
momentum that not even optimal choices of directions for the other two decays can bring 
Xi to rest (in the rest frame of qi)- For the decay which gives the maximum value of m qi u, 
the LSP will have a non-zero momentum in the same direction as the sparticle sent out 
from the 'hardest' decay. There are three separate dominance regions in mass space where 
this happens, one for each sparticle decay in the cascade. 

Consider the case where the second decay, that of is dominant. If Xi cannot be 
brought to rest in the rest frame of q~L, the maximum value is obtained for the optimised 
case where qi and If fly off in one direction and l n in the opposite, see Fig. ||. In the last 
step of the figure all momenta are boosted back to the rest frame of (Jl- If we are inside 
the dominance region, Xi ends up with a momentum pointing 'upwards' in the rest frame 
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Figure 8: One of three dominance regions for m 



max 

qtll 



of c[l, as in the figure. In this case the defining criterion of the region takes the form, 



IpJ > \p q + Pi f \ 



X2 y Ir "'■ 



my mr,o m^o 

lR X-1 X-2 



(A.4) 



where the three- momenta are taken in the rest frame of <jl. Note the form of the constraint 
in terms of ratios of the masses of nearest neighbours in the decay chain. The dominant 
ratio is the one which corresponds to the decay of %2- The expression for rn™f^ in this 
region is then found by simply combining the four-momenta of the three Standard Model 
particles for the given linear configuration. 

The defining inequality of the other two dominance regions and the resulting expres- 



sions for are found by similar approaches. The general result is again Eq. (4.4) of [12]. 

Note the systematic pattern of inequalities which define the dominance regions. Although 
in cases (1) and (3) the right-hand side of the inequalities can be simplified, this was not 
done in order to emphasise the structure of the constraints. 

The approach described above has some clear advantages over the more direct approach 
of method 1. Each region is 'understood' in terms of a dominant decay, here in the sense 
that the mass ratio is larger than the product of the other two, together with minimum 
momentum of Xi i n the rest frame of q^. No angle a is needed. Only cascade decays on 
a line need be considered, which allows the entire calculation to be performed on a few 
sheets of paper. 

Furthermore, and this is the most celebrated quality here, the approach is trivially 
extendable to an arbitrary number of particles. In particular, the solution for rn™5£ is 
found in the same way, and is given in Eq. ( |2.1| ). Notice the similarity with ni^f*. For a 
detailed calculation in the fully general case of a decay chain involving n sparticles, see [17]. 



B. Calculation of m^ low) 

The calculation of m ^m ow ) turns out to be quite involved. This is a consequence of the 
inherent complexity of calculating 'low'-endpoints, joined with the many possible decay 
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configurations available for the gluino decay chain. In order to demonstrate the peculiarities 
of 'low'-endpoints in a purer form, we first review the calculation of ^^^n ow ) > the simplest 
of the 'low'-endpoints. For the calculation of m q n q f^ ow ^ which follow thereafter many of the 
considerations are similar, only significantly more difficult to carry through. 

B.l Preparations: review JTi™fn w) 

Let us first recall that since there are two leptons, two invariant masses m„A can be con- 
structed, one of which will be higher than the other. These will be denoted the 'high' 
and 'low' distributions, each of which will have a maximum, denoted 'max'. Thus, the 
difficulty lies in identifying whether 'high' and 'low' correspond to the 'near' (l n ) and 'far' 
(Zf) leptons, or vice versa. The rn^n^^ and m qf i^ ow ^ distributions are constructed from 
the highest/lowest of m qf i n and m qi i { on an event by event basis. (We here assume that 
the correct jet has been selected.) Since the m qi i value which gives the absolute max- 
imum must necessarily be the higher of m qi i n and m qf i f for the given configuration, we 
simply have w™^ lg j^ = max(m™^ x , m™^ x ). For m^fn ow \ the situation is more complicated; 
we need to look for the maximum value of the lower of m Qi i n and m qt i r This conditional 
maximisation requires that both m qi i values are compared for the given configuration. 
Under no circumstance can the endpoint be higher than the lower of the two maxima, 

m m f? \ < minfrii" 1 ?* m™f x ). 
9ft (low) — v gftn ' q^l^ > 

Let us first assume a mass scenario where ml^f* < m ^ X i which directly corresponds 
to the condition 

rrir > m^om^o. (B.l) 



This is manifest in regions (1) and (2) of Eq. (4.5) of [12|. Now, let us further assume that 
we have a decay configuration where m q{ i n takes on its maximum value, m qf i n = rn^f*, 
so that (ft and l n are back-to-back in the rest frame of If it is possible to choose an 
orientation of p; f such that m qf i { > m qi i n , then m qi i a can be a 'low'-value, and we will have 
m q a [ i(iow) = m qltn' Specifically, this occurs when m qf i { takes on its maximum value, with 
p^. in the opposite direction to p q{ (and parallel to P/ n ), which corresponds to 

2mj >ml + ml , (B.2) 

Ir Xi X2 

i.e. region (1) of Eq. (4.5) of @. 

If this mass constraint is not satisfied, m q{ i n will be the 'high'-value for any configura- 
tion which has m q{ i n = m™f* and m qf i f will be the 'low' value. To find the configuration 
which gives m m f, x N in this case, consider the more general situation where the angle a 

gft (low) ' 

between p 9f and p; n (in the rest frame of X2) ^ s reduced. The maximum value of m q{ i f for 
this case is denoted m™^ x (a). [It will occur when p ?f and p l{ are back to back in the Ir 
rest frame since |pj | obviously is independent of the orientation in this particular frame.] 
If a is reduced (starting from n), m q{ i n (a) will decrease and rn^f^a) will increase. As a 
decreases, the 'high'-value decreases and the 'low'-value increases, and for some angle a eq 
they become equal: m qf i n (a eq ) = m™^ x (a eq ). If a is reduced below a eq , m qn i n will become 
the 'low'-value, and as it decreases with a the 'low'-value will decrease. Hence, a e q gives 
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the optimal configuration, ?™™fQ ow ) = m qf i n (a cq ) = m™ { i( eq y A simple calculation gives 
Eq. (4.8) of pi, and is valid for 

mlo + m 2 > 2m? , (B.3) 

•M A2 L R 



i.e. regions (2) and (3) of Eq. (4.5) of |L2|]. 



One might expect to also find a separate solution for m q{ i { < m qf i n and m qt i { = m™f*, 
in analogy with the first solution obtained above. The situation is however not symmetric 
with respect to the two leptons. Maximisation of m q{ i { necessarily fixes m qi i n at zero, so 
there is no third possibility, and once again m ™f(i ow ) = m |^/(eq)- The general solution for 
m <frZ l (iow) * s then given by two expressions, m™f* and m ^jfaq)i m the mass regions derived 

above, see Eq. (4.5) of (Below, we will see that a third possibility is relevant for 

max \ 
" l qql(l(my) 

B.2 m m ^ w introduction 

For a three-particle endpoint, the configurations which give the maximum values nt™ 1 ^ and 
m qqi*' as wen as the endpoint expressions themselves, will be mass dependent [contrary 
to the situation for m^f* and ni™^]. In particular, we see from Eq. (4.4) of [T^j and 



Eqs. (2.14) and ( |2.15| ) of this paper that there are four different cases (mass regions) for 
each. This increases the number of situations to consider, as well as the complexity for 
each. 

To calculate w^n ow ) we use the same strategy as above. We will first investigate the 
situation m qq i n < m qq i { with m qq i n = (Sect. |B.3| ) and find the appropriate conditions 

on the masses corresponding to w^n ow ) = m ^t*- Then we will investigate the 'opposite' 
situation, m qq i { < m qq i n with m qq i { = m,™| f x (Sect. |B.4| ), and find the mass conditions for 
the solution m-^Q ow ) = m q n q l i^- This differs from m^n ^ where no such solution was 
available. Finally an 'equal'-solution will be sought (Sect. |B.5| ). It will consist of a critical- 
point solution and boundary solutions (at the boundary of a two-dimensional parameter 
space). The calculation of the 'equal'-solution is cumbersome and preferably performed 
with the aid of a computer program. The resulting expressions for the full solution are not 
very complicated, but many. Care should be taken in the implementation. 

In the rest frame of the quark sector and the lepton sector can be described without 
reference to each other. This makes the rest frame particularly convenient to use. Unless 
stated otherwise, momenta and energies are given in this rest frame. The coordinate system 
will be chosen such that the combined momentum of the two quarks, p qq will point upwards 
(along +z). The x and y components are not required, allowing the four- vectors to be given 
solely by the energy and the momentum in the z-direction. 

13 o ™max ,™max 

±3 - 3 m qql(low) ~ m qql n 

We here investigate the situation m qq i n < m qq i { with m qq i n = rn 1 ^*. The possible config- 
urations with a maximised m qq i n are given in Table |l0[ The four rows correspond to the 
four different regions of mass-space appropriate to Eq. ( [2,14D (i.e. Eq. (4.4) of with 
the substitution of Eq. ( p,14[ )), as labeled by the first column. The second column shows 
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the corresponding 'region condition'. The third and fourth column show the directions of 
the quark and lepton momenta, from left to right as in the decay chain. The combined 
momentum of the quarks is always upwards (in the rest frame of x^)- The direction of 
p^ n is always downwards, while the direction of p^ f is upwards for configuration 1, down- 
wards for configuration 2. Notice that in region (4), configurations which correspond to a 
maximum of m qq in -w) do not have the quarks aligned in the x[> rest frame. There are four 
mass regions and two directions for p if , giving in total eight possibilities. Each must be 
considered. 





region condition 


config. 1 

<?n<?f Ink 


config. 2 
9n9f kik 


region (1) 


m iL m-o m- lR 


Tl IT 

superfluous 


TI II 

trivial 


region (2) 


m iL ^ m x° 2 m- g 
m x° 2 m l R mi} L 


IT IT 

superfluous 


IT II 

trivial 


region (3) 


m x% ^ m- g m iL 
m; m,g T m-n 


TT IT 

no solution 


TT II 

trivial 


region (4) 


otherwise 


\/ IT 

superfluous 


\t II 
trivial 



Table 10: Possible configurations for See the text for details. 



The comments 'trivial', 'no solution' and 'superfluous' give the conclusion of the in- 
vestigations. Some of these are seen right away: in region (3), configuration 1 the value of 
m qq [ f evidently vanishes since the three particles go in the same direction, hence m qq i { must 
always be smaller than m qq i n and there is 'no solution'. For configuration 2, in all four 
regions the condition m qq i n < m qq i { is simply \pi n \ < \Pi f \, which 'trivially' corresponds to 

m-o + m 2 - < 2m? . (B.4) 

The three regions of configuration 1 which are marked 'superfluous' give mass con- 
ditions which turn out to be subsets of the corresponding ('trivial') mass condition of 
configuration 2 [when the appropriate region condition is also imposed]. This means that 
if the ('superfluous') configuration 1 is possible, then also configuration 2 is possible. The 
conclusion of this subsection will therefore be that m ^fp ow ) = m ^t^ if anci only if the 
trivial condition flB.4| ) is satisfied. The rest of the subsection is dedicated to proving that 
the superfluous conditions are indeed contained in the trivial one. 

For each of the three regions we must first find the mass condition by requiring 
m qqi n < m qqi f - The proofs that these conditions are already contained in the trivial condi- 
tion [given that the appropriate region conditions are also imposed], are carried through by 

combining these with the 'anti-trivial' condition, m? n -|-m-n > 2m? , and see that this leads 

to ' xi xi Ir 

to contradictions. For regions (1) and (2) the proofs are fairly straightforward, but less so 
for region (4) due to the complicated nature of the region condition and the non-aligned 
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quarks. Below we first find the mass conditions of the three regions, then show the proof 
for one of them. The proofs of the two others are similar. 

Regions (1) and (2) 

In the rest frame of \% the relevant four-vectors of configuration 1 are: 

Pqq = (E qq , \p qq \), Pl n = (\PlJ, ~\PlJ), Plf = (\Pl f \, \Plf\) (B.5) 

where each entry is given in the form (E, p z ). The condition we need to satisfy is then: 

m qq l n < m qqk ^ (p qq +Pl n f < {p qq +Plff Pqq " Pl n < Pqq ' Pl { 

i E m + Ip m I)IpJ < ( E ii ~ |P OT l)|Pjfl ( B - 6 ) 
The lepton three-momenta have the magnitude: 

(m 2 -m 2 ) (m~ - m| )m-o 

|_ I X 2 l_R I I _ *1 *2 /-p 

|PJ " 2m,o ' |P ^- W (R7) 

x 2 2 H 

For region (1) the quark sector gives 

E m = |P<J + |Pg f l; |P OT I = |P fc l - IPj ( B - 8 ) 

+ ip«i = 2| Pfc i = !lJL -^-A ^ - Ip«i = 2 ipJ = - J ^ r ^ L - ( B - 9 ) 

X 2 

For region (2), the only difference is the expression for \p qq \, which differs by an overall sign 
(IPqijI = P<?J — Ip<?J)> an( ^ amounts to interchanging the two expressions in ( |B . 9| ) . If the 
expressions for regions (1) and (2) are inserted into ( |B.6| ), the following mass conditions 
are found: 



region (1): m~ (m~ — m~ L )(mio — m~ ) < m qL { m q L ~ m ^,o)( m f ~ m |o) (B.10) 

region (2): m~ L m~^(rn~ L - m|o)(m|o - n%~ ) < m|o(m| - m^ L )(m 2 R - m|o) (B.ll) 

If ( |B.10| ) or ( |B.11| ) is satisfied, together with the appropriate region condition (Table |10| 
second column), then we have m q ^n ow \ = m ^/*- 

Region (4) 

In region (4) the quarks are not sent out on a line (parallel/antiparallel). We must instead 
find p qq in the rest frame of X2 using the following equations: 

(m- g - m h ) 2 = p 2 qqln = p\ q + 2p qq ■ p ln = E 2 qq - p\ q + 2(E qq + |p M |)|pJ (B.12) 
m l = p lgx° 2 = P *<i + m l° 2 + 2p xl ' p ii = E li - p m + m \% + 2m x° 2 E ii ( R13 ) 
The first equation uses the expression for rn™*j£ found in Eq. ( |2.14 ). From these we find 

Pqq'- 



m g (m 2 +mf ) - 2mi mr m g {m 2 - mf ) 

£ _ _ X 2 lR X 2 L R jp I _ X 2 Ir ^\ 

qq 2m { omr ' qq 2m { omr 

X 2 tfl X 2 Ir 
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Eqq + |P 



m^o{m~ g - m~ lR ) 



qq I 



E, 



{m- g m~ lR 



mr 



qq 



J qq\ 



m x% 



Insertion of these expressions into ( |B.6[ ) gives the condition for region (4): 



,, - in? )(m- — m r ) < (msmr —m~ )(m~ .. ., 



m 2 ~ ) 



(B.15) 



(B.16) 



If this is satisfied together with the corresponding region condition in the second column 



of Table 10L then m™* ow) 



m 



qqln ' 



However, as stated earlier the mass conditions ( [B.IOD , ( B,ll| ) and (B.16) are 'superflu- 
ous'. Below this is proved for region (1). Prom the region condition and the 'anti-trivial' 
condition, we have the following inequalities: 



9 . m QL m X°2 _ _2 /_2 _2 ^ ^ _2 /_2 



> 



m 



X2 IR 



Z f Z Z \ ^ Z / 

<^ mr (m; — rrig ) > rrig (m; 



1L 
2 



Xl X2 *H v X2 Ir Ir X\ 

Starting from ( B.10| ) and using ( B.17| )-( B.18| ) we can then write 

2/2 2 \/ 2 2 \ ^ 2/2 2 \/ 2 2 \ 

m fc( m ?L " m x° 2 ^ m iR ~ m xV > m i R ( m 9 ~ m ?J( m x° " m i R > 

^ 2/2 2 \/ 2 2 \ 

> m q L ( m <iL ~ m i R )( m X l i ~ m I R > 
. 2/2 2\/2 2n 



(B.17) 
(B.18) 



mr > m^o 
Ir X2 



(B.19) 



which is clearly wrong. Hence, given the correct region condition, the mass condition (plc| ) 
is in contradiction to the anti-trivial condition, meaning that it must be contained in the 
trivial condition. Similar proofs apply for the two other superfluous cases. 



We therefore have the following conclusion, as already stated: if (B.4) is satisfied, then 



m qll(\aw) ~ m qqln ' 



.4 m , n \ — m„„, 

qql(low) qql{ 



Here we investigate the situation m qg / f < m gq / n with m gg i f = m^\ The possible con- 
figurations with a maximised m q9 / f are given in Table [ll] (in contrast to Table 10 where 
m qq i n was maximised) . As before we consider the rest frame of xH an d align the coordinate 
system such as to have p qq pointing upwards. The maximum value for m qq i { is found for 
the decay configuration of X2 an d Ir which maximises p^ f downwards (in the rest frame 
of ^2)- This is achieved if first Ir, then if, are emitted downwards. Since this necessarily 



fixes p; upwards, there is only one configuration in Table 11. The mass ratios in the 



quark sector are the same for m^f as for m™f*, but in the lepton sector the relevant ratio 



qql n ' 



is now m^/fm^m^/m^) = mj^/m^o, see ( p. 14 ) and (2.15), which changes the region 



conditions in column 2 as compared to those of Table 10 



Region (^J has no solution since m qq i n vanishes and therefore cannot be the 'high'-value. 
In the three other regions there are unique mass conditions for the solution m 



m 



qqh 



qql(\ow) 

which is what we seek here. The calculations follow the exact same path as in the 



previous subsection. 
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m qqh 


region condition 


config. 1 

?n?f kk 


region (1) 


m - 77}.- 77%7 

m g ^ il 1 r 
ma t m-o m-o 


T 1 T 1 

1 -4- 14- 

unique 


region (2) 


m Sr m ln ma 

1L \ l R 9 
m-o m-o mg L 


IT T I 

unique 


region (3) 


X{ q L X2 


TT TI 

no solution 


region (4) 


otherwise 


\t TI 

unique 



Table 11: Possible configurations for m™^. See the text for details. 



Regions (1) and (2) 

In the rest frame of the relevant four- vectors are [compare with ( |B.5| )1: 

p qq = {E qq , \p qq \), pi n = (|pj, [pj), p /f = (| P/f j, -|p Jf |) (B.20) 
The condition we need to satisfy is now: 

m qqlf < m qq l n (p qq +Pl f f < (p qq + Pl n f P qq ■ Pl { < P qq ■ Pl n 

O [E qq + |p M |)|p if | < (E qq - \p qq \)\p ln \ (B.21) 



The lepton momenta are still given by ( |B.7| ). The energy and momentum of the quarks 
are the same as in the previous subsection: Eq. ( |B.S| ) for region (1), the expressions inter- 
changed for region (2). 

When these expressions are inserted into ( |B.21 ), we find the following mass conditions: 



region (1): m|o(m| - m? i )(m? { - m|o) < m \ L m j R { m \ ~ m |»)( m |o ~ (B.22) 

region fH): m~ (m~ — m 2 - )(m~ — m? ) < m~ (m~ — m~ )(m? — "*? ) (B.23) 

If ( |B.22i ) or ( |B,23| ) is satisfied, together with the appropriate region condition (Table [TT| , 
second column), then we have m^jJ^^ = "i™^. 

Region 

In region the combination of 

(mg - m^m^o/m^) 2 = p^ f = p\ q + 2p gg • p Jf = E qq - p 2 qq + 2{E qq + |p w |)|p lf | (B.24) 
[equivalent to ( |B.12| )] with (B.13) results in: 



E qq - 7, 1 \Pq q \ - o [p. lb) 

2mr m^o Imr m^o 

Ir Xi Ir Xi 

(mgrnr — m^omjo) fmsmco — m^omr ) 

rp i |„ I _ v y X 2 Xi y T-! I 1 v y Xi X2 Ir' /r> o/;\ 

^gg + |P M | - — , E qq - \p qq \ - — (B.26) 

Xi Ir 
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Insertion of these expressions into ( |B.21| ) gives the condition for region (4): 
m 2 -o(m~ g m lR - m^om^m^ - m| ; ) < m^o (to^o - m^m~ lR )(m\v - mfj (B.27) 



max 



If this is satisfied together with the corresponding region condition, then ^^n ow ) = m qqi t ■ 
This result completes the special-case solutions of TO™|p ow y 

■R c m m£K _ max 
13 ° m qql(low) — m qql(cq) 

If none of the above mass conditions are fulfilled, the maximum of m qq n\ ov ,) will in general 
be reached in an 'equal-solution' decay configuration, where m qq i n = m gq i f . 

The relevant information on the two first steps of the cascade decay is contained in 
the value of \p qq \- (We stay in the rest frame of %2-) Its 

maximum value is attained if q n 

and (/f are sent off in the same direction, and is given by, 

2 2 
m~ — m~ 

IpSI— ^ (B.28) 
If sent off in opposite directions, \p qq \ is at its minimum, 

12 2 4 I 

m~mi — m~ 

Hi\ = -^4 (B.29) 

QL X2 

In the following we let jp J be a free variable and only later impose the physical constraint 



'qq\ 

|p w l g OpmIiIpmO- 

As before, we align the coordinate system to have p qq along +z. The first lepton, Z n , 
is emitted at an angle a relative to p qq . The direction of the second lepton, Zf , is chosen to 
maximise m qq i { for the given a, allowing the notation m^S^iot). We have this freedom since 
we are interested in the maximum, not just any rn qq i^ ec ^ value. At this stage the decay is 
parametrised in terms of two variables, \p qq \ and a. By requiring m qq i n (a) = TO^^a), the 
angle can be expressed in terms of \p qq \: 

_ m lJ m h + ^ ~ (\j m l + \P*l\ 2 ~ m X 2 )- m X% ( m I R ~ m lo)<P(\Pqq\) 



X 2 



(ml - m? ) (2m? - m 2 )\p 



(B.30) 



where </?(|p„J) is given by 



¥>(|P M I) = [ m lo|p OT | 2 + (m| ? + m|o - 2mfJ (m\ + m|o - 2m^m§ + |p gg | 2 )] 7 

(B.31) 

[Later we will need to ensure that | cos a oq [ < 1.] Insertion of this solution into Tn qq i n (a) [or 
m™^(a)] returns an expression for m qq i( eq ) in terms of \p qq \ only. Since we are interested 
in a maximum value we now search for the critical point of m qq i^ eq ^, and find an expression 
for the critical momentum, 



|p crit l 

Wqq I 



2m,- (m% + m? - 2m~ ) + (3m| - m\ ), /m? + 2m? 

l R Xi X2 'u' «ii Xi'V 9 i fl 



2 2 

mi n — m.„ 

Xi X2 



- 35 - 



a 2 2 2 

-4m7.m~ m r 
y X2 Ir 



X2 lR 



(B.32) 



Inserting this back into m, 



qql(eq) 



gives an expression for the critical point solution, 



m x% ^ m l R - m x? ) V 1 P ^ 1 2 + m S 

m k m l - K - m^M\p^\)] /(2ml ~ m I?) (B.33) 



All expressions found from the procedure above are formal. In order for them to be physical, 
the momentum and the angle a of the critical solution must lie in the allowed regions: 



Ip£*I e [IpJJMpBI 

_crit i 



COSOeqClp™]) G HM] 



(B.34) 
(B.35) 



For a given set of masses Eqs. (B.34) and ( B.35| ) must be tested numerically for the resulting 



-.crit I 
<?<? I 



If the critical solution is not physical, the maximum must lie on the boundary of 
the two-dimensional domain defined by Eqs. (B.34) and ( B.35[ ). For cosa eq = ±1 the 
appropriate \p qq \ can be found by redoing the procedure leading to ( B.33 ) but with a fixed 
at or 7r. This gives the following expression for |p„J: 



I cos a=±l I 
Wqq 



±(m 2 m 2 - mf ) 

Xi Kl LR 



im^omf (m| - mf )(mf - m| ) L 

A2 l R X2 Ir Ir Xi 



+ 



\l {rn 2 - m% - mf ) 2 + 4m~mj (ml - mf )(m~ - ml 

V Xi X2 Ir y Ir X2 Ir Ir Xi 



/n 2 2 2 2 4 \ 

(2m~ m T —m~ m~ —m T 
y X" Ir Xi x" 

(B.36) 



The invariant mass can then be found straightforward by taking Tn 2 ql ^ = (p qq + pi n ) 2 and 
using ( B.13 ). This gives 

,cos «=±] \ 2 



[m 



qql(eq) 



K' *2 ~ X 2 
2 m„cosq=±1 



R 

l m x\ 



(B.37) 



The boundary solutions at maximum and minimum |p M | are most easily found by 
inserting the boundary values ( |B,2g| )-( [B.29 ) into the general solution (B.33). This gives 



Ki (eq) ) = (™1 - - m h)/( 2m I R - m h) (B.38) 

Ki(c q )) 2 = [ - ( 3m i R - m l°)( m l m l§ + m fJ + 2m lA m §( 2m i R - m l°) + m l° m2 J 

+2m] L {m 2 lR - m|o)^(| P ^|)]/[2m? L (2m? B - m| ? )] (B.39) 

All expressions ( |B.37| )~( |B.39| ) are formal. For a given set of masses one must explicitly 
(numerically) impose the one of Eqs. ( B.34 ) and ( B.35| ) which is not satisfied by construc- 
tion. 
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B.6 General solution 



Putting all this together, the fully general solution for m 



rn 



qgi(low) 



m 
rn 
rn 



gqi(low) 
2 j_ ™2 ^ 2 



is given by 



max 



qql 



max 
qql(eq) 



for mtn + m~ n < 2m~ 

x" X2 i R 

for (cutl) 
otherwise 



where 
(cutl) 



rn 



!J_ > ^HlJ^R_ A 



m 



X2 



m- 



m- 



qi 



m 



Ir 



x? 



1-lL 



X-2 Xl 



m? m~ [ml 



2 

X2 



2 2 

mf n — m~ 
Xi Ir 



< 1 



V 
V 



m 



<1L 



m 



X-2 



> 



9 

X2 



2 



2 

x? 



to- ( to- 



to 



Hi 



m 



Ir 



< 1 



< 



rn 



m q L m iR A m 9L 



< 



X-2 Xl 



A.2 



A 



S < m 9 



Xl 



mq L m^o 



A 



2 

x? 



»i? Xi v y Xi X2 'h 



x' 2 ' 



(B.40) 



< 1 



(B.41) 



is constructed from (B.22), flB.23j ) and ( |B.27| ) together with the appropriate region condi- 
tions. Expressions for to™|* and to™^* are given in Eqs. (2.14) and ( 2.15) ) respectively, and 



to. 



°qql(eq) 

which uses (|R3|) and (|Bl^)-( lR39|) . 



l,,aX I m g5eq)» m ^(eq) ± ^ m JgKeq)' 77t ?^(eq) 



(B.42) 
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